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ABSTRACT 


A  frequency-domain  method  of  predicting  the  incidence  of  personnel 
or  equipment  sliding  is  presented  which  includes  the  (linearized)  forces 
due  to  roll,  pitch,  longitudinal,  lateral  and  vertical  accelerations,  and  also 
the  effects  of  non-zero  mean  heel.  The  predictions  of  the  method  are 
compared  with  the  results  of  observations  of  at-sea  sliding  incidents,  and 
good  correspondence  is  obtained. 

A  simple  approximation  for  the  duration  of  threshold  exceedance  prob¬ 
ability  density  is  described,  and  used  to  develop  a  method  for  predicting 
the  incidence  of  slides  greater  than  a  certain  severity.  The  method  applies 
for  arbitrary  lateral  and  vertical  acceleration,  but  only  to  the  case  in  which 
the  longitudinal  acceleration  is  negligible.  An  example  of  the  application 
of  the  method  is  presented. 

It  is  argued  that  ship  motions  such  as  roll,  pitch,  lateral  and  longitu- 
din.\l  accelerations  are  the  wrong  physical  parameters  for  expressing  sea¬ 
keeping  criteria.  These  criteria  become  simpler  when  expressed  in  terms 
of  the  incidence  .->f  degrading  events  such  as  helicopters  sliding,  and  are 
equally  applicable  to  all  vessel  types. 

ADMINISTRATIVE  INFORMATION 

This  work  was  done  while  the  first  author  was  on  exchange  at  the  David  Taylor  Re¬ 
search  Center  from  Defence  Research  Establishment  Atlantic.  It  wew  sponsored  by  the 
Chief  of  Naval  Research,  Office  of  Naval  Technology,  Code  0NT211,  under  6.2  Surface 
Ship  Technology  Program  (NDIA),  Program  Element  62121N,  Advanced  Hull  Project 
RH21S23,  Task  5,  Ship  Motion  Control.  The  work  was  performed  during  FY1991  under 
work  unit  number  1-1606-122.  The  DN  number  is  DN178067. 

INTRODUCTION 

Predictions  of  the  operational  performance  of  ships  at  sea  require  three  elements;  a 
statistical  description  of  the  wave  and  wind  environment  in  which  the  ship  will  oper¬ 
ate;  computer  programs  for  determining  the  ship  response;  and  seakeeping  criteria  to 
indicate  the  effects  of  the  motions  on  operational  performance.  Hindcast  data  bases  are 
now  available  for  many  operational  area.s  (see  e.g.  Reference  1),  and  strip  theory  com¬ 
puter  programs  such  as  SMP^"*  are  in  wide-spread  use.  There  has  been  comparatively 
less  effort  devoted  to  developing  seakeeping  criteria,  and  the  accuracy  of  operability 
predictions  is  now  limited  by  the  accuracy  of  available  criteria. 

A  summary  of  criteria  in  use  in  1989  was  presented  in  Reference  4.  As  an  example, 
criteria  for  personnel  consist  of  separate  limits  for  each  of  the  following  four  parameters: 
roll,  pitch,  and  lateral  and  vertical  acceleration  at  the  pilot  house.  There  are  some  prob¬ 
lems  associated  with  this  approach.  First,  operations  are  not  limited  by  any  one  of  these 
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parameters  individually,  but  rather  by  some  combination  of  all  four.  A  value  of  roll 
which  will  cause  extreme  problems  in  the  presence  of  significant  vertical  accelerations, 
will  be  perfectly  acceptable  when  vertical  motions  are  small.  Second,  the  criteria  are 
primarily  based  on  operational  experience  with  frigates  and  destroyers.  While  the  lim¬ 
iting  values  of  these  parameters  found  do  reflect  the  level  of  motions  at  which  significant 
personnel  problems  occur  on  these  vessels,  they  may  result  in  misleading  predictions 
of  the  operability  of  SWATH  vessels  or  even  of  monohulls  of  significantly  different  size. 
For  example,  the  pitch  criterion  of  LS""  RMS  w^hich  is  based  on  operational  experience 
on  vessels  with  length  to  beam  ratios  of  about  9  with  most  frequently  occurring  pitch 
periods  of  around  6  seconds  might  be  conservative  if  applied  to  a  SWATH  vessel  with 
a  length  to  beam  ratio  of  2  and  longer  pitch  periods. 

In  fact,  roll,  pitch,  vertical  and  lateral  accelerations  are  the  wrong  physical  param- 
eters  for  expressing  personnel  criteria.  The  limiting  values  of  these  parameters  found 
through  operational  experience  are  representative  of  the  motion  levels  at  w^hich  certain 
degrading  effects  start  to  become  important.  For  personnel  operating  on  deck,  the  de¬ 
grading  effects  are  the  onset  of  loss-of- balance  events  due  to  tipping  oi  sliding.  These 
events  were  called  motion-induced  interruptions  (Mil)  in  Reference  5.  Reference  6  in¬ 
troduced  a  frequency-domain  method  for  predicting  the  incidence  of  motion-induced 
interruptions  that  included  the  effects  of  both  lateral  and  vertical  acceleration,  and 
proposed  that  a  rational  seakeeping  criterion  for  deck  operations  involving  personnel 
could  be  developed  in  terms  of  the  incidence  of  MIFs  per  unit  time. 

In  future,  we  propose  that  all  seakeeping  criteria  should  be  presented  in  terms  of 
limits  on  the  incidence  of  whatever  degrading  event  is  determined  to  limit  the  operation 
of  interest  rather  than  the  underlying  ship  motions.  As  an  exan^ple,  we  consider  an 
unsecured  helicopter  on  deck,  as  addressed  in  Reference  7.  In  this  case,  the  degrading 
event  of  most  interest  is  the  onset  of  helicopter  sliding.  A  criterion  in  terms  of  the 
maximum  acceptable  number  of  slides  per  hour  is  physically  easy  to  understand,  and 
obviously  applies  equally  to  SWATH  or  monohull  vessels.  The  alternative  approach  of 
trying  to  determine  separate  limits  for  the  contributing  motions  of  roll,  pilch,  vertical, 
lateral,  and  longitudinal  accelerations  as  well  as  relative  wind  will  not  lead  to  a  unique 
answer  because  of  their  interrelation;  moreover,  different  values  w^ill  be  obtained  de¬ 
pending  on  whether  the  deck  is  slippery  or  not.  If,  on  the  other  hand,  the  criterion  is 
expressed  in  ieims  of  a  limit  on  the  allowable  number  of  slides  per  unit  time,  the  effect 
of  the  coefficient  of  friction  is  automatically  taken  into  account.  It  should  be  noted  that 
the  philosophy  outlined  above  is  already  used  routinely  for  slamming  and  deck  wetness 
criteria  which  are  usually  expressed  2is  limits  on  the  number  of  slams  or  wetnesses  per 
hour. 

In  order  to  implement  our  approach  to  developing  seakeeping  criteria,  two  steps  are 
required.  First  methods  must  be  developed  to  predict  the  incidence  and  severity  of 
degrading  events  such  as  sliding  of  helicopters  or  MIFs  of  f>ersonnel  on  deck  In  order 
to  lead  to  practical  tools  for  ship  design,  these  methods  must  w'ork  in  the  frequency 
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domain.  Second,  operational  data  are  required  to  determine  acceptable  limits  for  the 
incidence  of  occurrence  of  the  degrading  events. 

In  this  report,  which  is  an  extension  of  Reference  6,  we  present  a  frequency-domain 
method  of  predicting  the  incidence  and  severity  of  personnel  or  equipment  sliding  which 
includes  the  (linearized)  forces  due  to  roll,  pitch,  longitudinal,  lateral  and  vertical  ac¬ 
celerations,  and  also  the  effects  of  non-zero  mean  heel  angle.  An  analysis  of  personnel 
tipping  which  includes  the  above  forces  and  also  steady  and  unsteady  wind  effects  will 
be  published  later*. 

PREDICTION  OF  SLIDING  INCIDENCE 
WITH  ZERO  LONGITUDINAL  ACCELERATION 

In  this  section,  we  review  the  frequency-domain  method  for  predicting  sliding  in¬ 
cidence  in  the  absence  of  longitudinal  acceleration  which  was  presented  in  Reference 
6.  Longitudinal  acceleration  is  usually  negligible  for  large  monohulls,  but  need  not  be 
for  small  monohulls  or  different  platform  types.  A  method  for  predicting  sliding  inci¬ 
dence  in  the  presence  of  combined  longitudinal,  lateral,  and  vertical  acceleration  will 
be  presented  in  a  later  section. 

Two  coordinate  systems  wdll  be  used  for  the  sliding  calculations.  The  inertial  system 
is  illustrated  in  Fig.  1.  Its  origin  is  located  at  the  mean  position  of  the  ship  and  it 
translates  with  the  mean  velocity  of  the  ship,  maintaining  a  fixed  orientation  with 
respect  to  the  free  surface.  The  motions  surge,  sway,  heave,  roll,  pitch,  yaw  will  be 
denoted  by  r/,,  i  =  1,  •  •  •  ,6. 

The  displacement  D  =  {Di,  D2^  at  a  point  P  =  (JV,  K,  Z)  on  the  ship  is  given 
by 

D  =  if+ixP  (1) 

where  ff  =  (771,772,7/3)  and  ^  =  (^1,^2, ^3).  For  example,  the  displacement  in  the  vertical 
direction  is  given  by 

r>3  =  773  +  (2) 

and  includes  contributions  from  heave,  roll,  and  pitch.  The  velocities  and  accelerations 
at  P  are  obtained  by  differentiating  Eq.  1  with  respect  to  time. 

3  =  (3) 

3  —  rf  +  [  X  P  (4) 

*Graham.  R.,  Baitis,  A.E.,  and  Meyers,  W.G..  “On  the  Development  of  Operability  Criteria’*,  to 
appear  in  Naval  Engineers  Journal,  ASNE  Day.  1992 
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An  object  of  mass  m  fixed  to  the  ship  with  center  of  gravity  located  at  the  point  P 
experiences  an  inertial  force,  which  in  the  earth  system  is  given  by 

Fc  =  —  D2  +5)  (5) 

where  g  is  the  acceleration  due  to  gravity.  The  object  could  be  an  unsecured  helicopter 
on  deck,  a  crew  member,  or  anything  else.  For  the  development  of  seakeeping  criteria, 
the  forces  perpendicular  and  parallel  to  the  deck  are  of  more  relevance.  These  are  the 
forces  in  the  ship  reference  system,  F^. 

Fig.  1  also  illustrates  the  ship  coordinate  system,  which  is  now  fixed  to  the  ship. 
Performing  the  coordinate  transformation  from  the  earth  to  ship  reference  system,  and 
retaining  only  the  linear  terms  we  obtain 

F,  =  m[-Di  +  gT)s,  -£>2  -  gT]i,  -Z?3  -  g)  =  rn{Fiong.Fi^t^  F^)  (6) 

The  quantities  Fiong  and  Fiat  are  the  (linearized)  longitudinal  and  lateral  forces  per 
unit  mass,  and  following  Reference  5  will  be  called  the  longitudinal  and  lateral  force 
estimators.  F^  is  the  (linearized)  vertical  force  per  unit  mass. 

Suppose  that  the  object  under  consideration  has  a  static  coefficient  of  friction  given 
by  /i,.  Under  the  assumption  that  Fiong  ^  0,  slides  can  only  occur  in  the  port  or 
starboard  directions.  A  slide  to  port  will  occur  whenever 

^nFiat  >  -mfj^F^,  (7) 

or 

-D2  -  gm  -  Mi^3  >  ikg  (8) 

On  the  other  hand,  the  object  w'ill  slide  to  starboard  if 

mFjaf  <  (9) 

or 

Dz  +  5774  -  M  (10) 

The  quantities  on  the  left  hand  sides  of  Eqs.  8  and  10  were  called  generalized  lateral 
force  estimators  in  Reference  6.  In  the  present  work,  these  quantities  will  be  called 
sliding  estimator  functions,  and  denoted  by  5p  and  5,,  respectively. 

Wind  forces  on  the  ship  or  asymmetric  loading  conditions  may  result  in  the  presence 
of  a  steady  heel  angle,  (positive  starboard  side  down).  In  this  case,  the  right  hand 
sides  of  Eqs.  8  and  10  should  be  replaced  by  jj^g  +  4>rj^g  and  745  —  0mg,  respectively. 
The  net  result  is  to  make  sliding  more  likely  to  occur  in  the  downhill  direction. 

In  Reference  6.  a  frequency-domain  method  w’as  derived  for  predicting  the  number  of 
slides  per  unit  time,  under  the  assumption  that  the  sliding  estimator  functions  followed 
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the  Rayleigh  distribution.  In  fact,  the  result  holds  provided  that  Sp(t)  and  Sj(0  are 
zero-mean  Gaussian  processes,  and  is  not  restricted  to  narrow-band  processes.  Let 
Ss,{<*>e)  denote  the  (one-sided)  power  spectral  density  of  S,{t)  as  a  function  of  encounter 
frequency  w*,  and  denote  by  7n„  the  spectral  moments  of  S5,(u’r). 


m, 


=  [  w"5s,(u,',)cL.'e 

Jo 


(11) 


Rice®  (see  also  Reference  9)  showed  that  the  number  of  upcrossings  (or  downcrossings) 
per  unit  time  of  a  threshold  at  level  a  is  given  by 


1  I 

=  i^J—exp  (  - 
Ztt  y  ttiq 


2mc 


(12) 


The  quantity  M,  is  the  number  of  slides  per  unit  time  in  the  starboard  direction.  In 
the  absence  of  steady  heel,  the  threshold  a  is  at  iJ,g.  In  the  general  Ccise,  the  threshold 
a  is  at  The  number  of  slides  in  the  port  direction  Mp  can  be  estimated  in 

a  similar  m2inner. 

The  onset  of  widespread  sliding  of  objects  on  board  ship  is  indicative  of  a  level  of 
motion  which  makes  all  shipboard  operations  difficult,  and  necessitates  tying  down  all 
moveable  objects.  It  is  suggested  that  a  minimal  seakeeping  criterion  for  almost  any 
operation  is  that  the  probability  of  sliding  remain  r  cceptably  low. 

The  probability  of  sliding  depends  on  the  precise  location  on  the  ship;  hence,  in 
specifying  a  sliding  criterion,  it  is  also  necessary  to  specify  a  location  or  series  of  loca¬ 
tions  of  interest.  For  a  transit  mission,  for  example,  the  locations  of  interest  could  be 
the  bridge,  the  engine  room,  etc.,  and  the  sliding  criterion  would  have  to  be  evaluated 
simultaneously  at  the  various  locations.  The  incidence  and  severity  of  sliding  which  is  ^ 
considered  tolerable  at  each  location  will  ultimately  have  to  depend  on  the  tasks  be¬ 
ing  performed  and  the  degree  of  degradation  in  task  performance  which  sliding  causes. 
Normally,  this  degradation  will  result  from  the  reduced  pace  of  task  performance,  but 
in  some  caises,  e.g.  aircraft  sliding  on  deck,  sliding  may  result  in  injuries  to  personnel 
or  damage  to  equipment. 


COMPARISON  OF  MEASURED  AND 
PREDICTED  SLIDING  INCIDENCE 

In  order  to  validate  the  method  presented  above  for  predicting  sliding  incidence,  the 
results  of  a  chair-sliding  experiment  conducted  in  1989  on  board  the  USCG  Morganthau 
were  analysed.  Unlike  helicopter  sliding,  chair  sliding  is  more  likely  to  be  a  nuisance 
than  to  result  in  serious  degradation  of  operational  performance;  however,  the  chair 
sliding  experiment  described  here  provides  a  valid  check  of  the  method,  and  had  the 
advantage  of  being  easy  to  conduct  without  appreciable  danger  to  the  crew  or  equip¬ 
ment.  Moreover,  the  lower  coefficients  of  friction  on  the  inside  of  the  ship  allow  this 
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type  of  experiment  to  be  conducted  under  much  less  severe  conditions  than  would  be 
necessary  to  cause  an  appreciable  incidence  of  helicopter  sliding  events. 

Basic  ship  data  for  the  Morganthau  are  given  in  Table  1.  In  this  experiment,  a 
subject  sitting  in  a  standard  government  issue  chair  with  metal  feet  was  located  in  a 
compartment  with  a  linoleum  floor  which  had  been  pitted  and  scarred  from  heavy  use. 
The  subject  was  facing  to  port  and  was  located  in  front  of  a  computer  terminal.  At 
the  onset  of  the  chair  sliding,  the  subject  was  instructed  to  press  a  function  key  which 
time-stamped  the  data.  A  total  of  25  slides  were  recorded  during  the  course  of  the 
experiment,  which  lasted  20  minutes.  The  ship  wais  heeled  to  starboard  by  3.3  degrees, 
and  as  a  consequence,  all  of  the  slides  occurred  in  this  direction. 

Table  1.  Basic  Ship  Data  for  the  USCG  Morganthau 


Displacement 

3055. 

tonnes 

LBP 

106.7 

metres 

Beam 

12.80 

metres 

Draft 

4.40 

metres 

The  chair  was  located  about  1.5  metres  to  port  of  the  centreline  and  2.0  metres  aft 
of  the  ship’s  centre  of  gravity.  The  ship  motion  instrumentation,  which  was  installed  for 
an  unrelated  experiment,  w'as  located  near  the  flight  deck  at  a  centreline  position  26.9 
metres  aft  of  the  CG.  The  flight  deck  lateral  and  vertical  accelerations  were  corrected 
to  the  position  of  the  chair,  by  using  the  measured  roll,  pitch  and  yaw  signals  to  take 
into  account  the  difference  in  moment  arms.  The  longitudinal  acceleration  at  the  chair 
was  estimated  using  the  recorded  pitch  signal  and  the  distance  of  the  chair  from  the 
w'aterline  and  it  was  found  to  be  negligible. 

From  Eqs.  7  and  9,  slides  should  occur  whenever  Fut  >  —f^Fy  or  Fiat  < 
Unfortunately,  the  coefficient  of  friction  of  the  chair  on  the  floor  was  not  measured 
directly  during  the  experiment,  but  estimates  were  obtained  by  computing  the  values 
of  the  ratio  |F/a<|/|Ft,|  for  the  observed  sliding  incidents.  The  values  obtained  are  shown 
in  Table  2. 

In  practice,  only  17  of  the  25  sliding  incidents  were  included  in  Table  2,  since  it 
appeared  that  the  subject  was  slow  in  time-stamping  the  remaining  slides.  Motion 
data  was  recorded  at  three  samples  per  second,  and  so,  the  uncertainty  in  the  starting 
time  of  the  slides  was  probably  about  0.5  seconds.  The  values  shown  in  Table  2  result  in 
a  mean  coefficient  of  friction  of  0.189  with  a  standard  deviation  of  0.024.  The  variation 
in  the  ^  values  found  is  consistent  with  the  poor  condition  of  the  floor.  It  should 
be  noted  that  even  under  laboratory  conditions,  significant  variation  in  coefficients  of 
friction  are  found  when  measurements  are  repeated.  (See  for  example  Reference  10.) 

Figure  2  illustrates  a  sample  plot  of  Fiat.  and  -^F^j  versus  time.  From  Eqs.  7 
and  9,  slides  should  occur  whenever  Fiat  crosses  eitl.cr  of  the  other  two  curves.  The  time 
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Table  2.  Computed  Coefficients  of  Friction 


Time  of  Event 

/4 

Time  of  Event 

34.9 

.184 

967.3 

.191 

54.6 

.182 

1038.9 

.219 

75.5 

.205 

1117.2 

.207 

302.6 

.174 

1125.8 

.175 

500.1 

.174 

1141.8 

.176 

559.7 

.166 

1152.5 

.211 

776.8 

.168 

1176.1 

.211 

803.8 

.162 

1186.J 

.245 

827.1 

.161 

stamps  corresponding  to  the  recorded  slides  are  shown  as  vertical  lines  on  the  figure. 

The  same  information  can  be  presented  in  terms  of  the  starboard  sliding  estimator 
function,  which,  from  Eq.  10,  is  given  by  5,  =  —Fiat  —  /4-D3.  In  this  context,  Fut  is 
simply  minus  the  lateral  acceleration  at  the  chair  since  the  acceleration  was  measured 
in  the  ship  reference  system  and  therefore  already  includes  the  component  of  gravity 
parallel  to  the  deck  due  to  the  roll  angle.  The  mean  offset  in  Fiat  due  to  the  nonzero  heel 
angle  should  be  removed,  so  that  the  upcrossings  of  5,  can  be  estimated  using  Eq.  12, 
which  is  only  applicable  to  zero-mean  processes.  Slides  are  predicted  to  occur  whenever 
S,  exceeds  the  threshold  —  g4>m-  Figure  3  shows  a  plot  of  the  starboard  sliding 
estimator  function  for  th  i  same  time  interval  as  shown  in  Fig.  2.  The  threshold  f^g—g<Pm 
is  also  shown  in  the  figure.  Although  Figs.  2  and  3  contain  the  same  information,  the 
presentation  of  Fig.  3  results  in  a  constant  threshold,  and  allows  the  sliding  incidence 
to  be  estimated  in  the  frequency  domain. 

Figures  4  and  5  show  plots  of  the  sliding  estimator  function  for  two  other  time 
slices.  Good  correspondence  between  observed  sliding  incidents  and  threshold  crossings 
is  obtained  in  Fig.  4,  except  for  the  fifth  slide,  which  appears  to  have  been  noted 
late,  and  the  eighth  slide,  which  seems  to  occur  at  an  anomalously  low  S,  value.  The 
experimenter  noted  that  the  slides  at  the  end  of  this  time  slice  were  rather  large.  This 
is  in  qualitative  agreement  with  the  large  threshold  exceedances  which  occured  during 
the  sixth  and  seventh  slides. 

On  the  other  hand,  the  correspondence  between  threshold  crossings  and  observed 
slides  shown  in  Fig.  5  is  poor.  In  this  case,  the  two  events  shown  do  not  correspond 
to  threshold  crossings,  while  a  subsequent  crossing  did  not  result  in  an  event.  Figure 
6  presents  the  same  time  slice,  with  the  addition  of  the  starboard  sliding  estimator 
functions  and  thresholds  corresponding  to  coefficients  of  friction  of  plus  or  minus  one 
standard  deviation  about  the  mean.  The  effect  of  the  change  in  coefficient  of  friction 
on  the  sliding  estimators  is  small;  however,  the  effect  on  the  threshold  is  significant. 


The  two  events  shown  now  correspond  to  crossings  of  the  lower  threshold,  while  the 
‘non-event’  at  8S2  seconds  did  not  exceed  the  upper  threshold.  Effectively,  the  variation 
in  the  coefficient  of  friction  results  in  uncertainty  as  to  whether  a  particular  crossing 
will  result  in  an  event.  For  the  purposes  of  validating  the  method  presented  herein,  it 
is  most  important  that  the  average  number  of  events  be  predicted  correctly. 

Overall,  the  correspondence  between  threshold  crossings  and  observed  sliding  events 
was  as  follows.  Sixteen  of  the  observed  slides  corresponded  to  threshold  crossings, 
whereas  nine  events  did  not.  On  the  other  hand,  there  were  seven  upcrossings  which 
did  not  result  in  events  for  a  total  of  twenty-three  upcrossings.  This  corresponds  well 
with  the  observed  number  of  twenty-five  slides.  The  number  of  upcrossings  was  also 
estimated  using  Eq.  12  and  the  computed  standard  deviations  of  the  sliding  estimator 
function  and  its  time  derivative.  This  resulted  in  an  estimate  of  24.5  upcrossings  for  the 
twenty  minute  period  which  is  also  in  close  correspondence  w'ith  the  observed  number 
of  slides. 

It  should  be  noted  that  real  data  may  show  significant  variability  in  the  number 
of  upcrossings  in  adjacent  time  slices.  It  must  be  emphasized  that  the  method  pre 
sented  here  is  statistical,  and  only  determines  the  average  number  of  upcrossings  over 
a  sufficiently  long  period. 

PREDICTION  OF  SLIDING  SEVERITY  WITH 
ZERO  LONGITUDINAL  ACCELERATION 

For  some  applications,  notably  fixed-wing  operations  on  aircraft  carriers,  operators 
routinely  tolerate  numerous  small  sliding  incidents,  and  only  curtail  operations  when 
severe  slides  start  to  occur.  This  suggests  that  it  may  be  of  interest  to  predict  not 
only  the  onset  of  sliding,  but  also  to  predict  the  severity  of  sliding,  as  measured  by  the 
distance  that  the  object  slides  before  coming  to  rest.  The  purpose  of  this  section  is 
to  present  an  approximate  frequency- domain  method  for  estimating  the  probability  of 
occurence  of  sliding  incidents  of  arbitrary  severity.  The  method  only  applies  to  the  case 
in  which  the  longitudinal  acceleration  can  be  neglected. 

REGULAR  WAVE  CASE 

In  this  section,  we  treat  the  special  case  in  which  the  starboard  (or  port)  sliding 
estimator  function  is  sinusoidal,  and  we  first  assume  that  the  dynamic  coefficient  of 
friction,  /^,  and  the  static  coefficient  of  friction,  /i,.  are  equal.  In  this  case 

5,(t)  =  5,0  sinu;t  (13) 

and  the  situation  is  as  illustrated  in  Fig.  7.  VVe  denote  by  f*  the  time  at  which 

5,osinu;t‘  =  a  (14) 

where  a  =  -  0mg-  Choosing  t  =  t‘  &s  the  new  time  origin,  the  mass  experiences  a 

force  given  by 
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m(Srt  sin  uf(< +  #*)  — a^)  for  <  > 

Denoting  the  distance  in  the  starboard  direction  by  x,,  we  have 


(15) 


mx,{t)  =  m{Sto  sinu;(t  +  t")  —  ag) 


The  velocity  at  time  t  is  given  by 


i.(0 


-  cosu;<)  +  —  sinwt  —  at 
u 


(16) 


(17) 


for  times  t  >  t’  such  that  i,(0  is  greater  than  zero.  Once  the  mass  comes  to  rest, 
the  frictional  force  keeps  it  at  rest  until  the  next  threshold  exceedance.  Eq.  17  can  be 
solved  iteratively  to  determine  the  time  at  which  the  m2iss  stops.  Choosing  our  origin 
at  the  initial  rest  position  of  the  object,  the  distance  at  time  t  is  given  by 


i,(t)  = 


/  s\nu>t\  a  ,, 

V - Z~J  ^  ”  cosut)  - 


(18) 


The  distance  travelled  can  now  be  determined  by  substituting  t  =  <,  in  the  previous 
equation. 

In  the  case  Mi  9^  '•h®  force  per  unit  mass  on  the  object  after  the  onset  of  sliding 

in  the  starboard  direction  becomes 


- Fict  -  IJdDs  -  fJd9  +  <Pmg  (19) 

when  the  sliding  estimator  function  exceeds  the  threshold  (/4  —  4>m)9-  Equation  19  will 
be  approximated  by 

-Fut  -  -  (m<  -  <t>m)9  =  ■?»(0  -  (Mi  -  <i>m)g  (20) 


This  approximation  is  required  because  of  the  method  used  to  estimate  the  amplitude 
of  the  sliding  estimator  functions  in  the  frequency  domain  which  will  be  discussed  in 
the  next  section.  It  can  be  expected  to  be  good  whenever  the  vertical  acceleration  is 
small  compared  to  9,  or  when  /4  —  Mi  is  small. 

With  the  approximation  of  Eq  20,  and  a  sinusoidal  sliding  estimator  function,  the 
velocity  and  distance  travelled  in  the  starboard  direction  become 


i.(0 


( 1  -  cos u;<)  +  —  sin ut  -  {m  -  4>-m )9t 
u> 


(21) 


i,(')  = 


\/(5'io  -  /  sinwA  a  1  ,  .  ^  .2 

- -  t - +  — (1  -  cosu;^)  -  -(Mi  -  4>m)9i 

u;  V  uJ  /  I 


(22) 
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While  the  incidence  of  sliding  depends  only  on  the  sliding  distance  depends 
mainly  on  which  determines  the  magnitude  of  the  frictional  force  once  motion  has 
started.  As  an  example,  for  an  angular  frequency  of  0.628  (corresponding  to  a  ten- 
second  period  for  the  sliding  estimator  function),  an  amplitude  of  0.5  g's,  and  zero  heel 
angle,  Eq.  22  predicts  a  sliding  distance  of  0.57  metres  for  /j,  =  0.45  =  /<d,  but  a  sliding 
distance  of  4.36  metres  for  =  0.45  and  =  0.35.  Another  consequence  of  Eq.  22  is 
that  if  there  is  a  large  difference  between  /j,  and  there  are  no  insignificant  slides:  all 
slides  are  serious  ones. 


IRREGULAR  WAVE  CASE 

A  time  history  of  the  starboard  sliding  estimator  function  derived  from  the  chair 
sliding  experiment  is  shown  in  Fig.  8.  We  will  assume  that  in  the  vicinity  of  a  threshold 
exceedance  the  sliding  estimator  function  can  be  adequately  approximated  by  a  sine 
wave.  It  remains  to  develop  frequency-domain  methods  for  estimating  tne  amplitude 
and  angular  frequency  of  typical  threshold  exceedances. 

We  will  make  use  of  a  high-threshold  approximation  for  the  duration  of  exceedance 
probability  density  function  due  to  Rice”.  In  Appendix  A,  the  applicability  of  this 
approximation  to  the  wave  elevation  process  is  investigated,  and  is  shown  to  give  good 
results  for  thresholds  greater  than  or  equal  to  one  standard  deviation.  This  condition 
is  not  unduly  restrictive  for  present  purposes,  since  if  the  standard  deviation  of  the 
sliding  estimator  function  is  greater  than  or  equal  to  the  sliding  threshold,  it  is  clear 
that  severe  sliding  is  so  common  that  operations  would  be  impossible. 

As  before,  we  assume  that  the  sliding  estimator  function  5,(<)  is  a  zero- mean,  Gaus¬ 
sian,  stochastic  process  (not  necessarily  narrow-banded),  and  suppose  that  the  threshold 
level  is  at  a.  We  define  the  nondimensional  threshold  by  r  =  a/ a,,  where  a,  =  y/m^  is 
the  standard  deviation  of  S,(t). 

W'e  will  denote  by  r  the  duration  of  a  given  threshold  exceedance,  and  the  dura¬ 
tion  probability  density  by  p,{t).  From  Eq.  61  of  Appendix  A,  we  have  the  following 
approximation  for  the  duration  probability  density  p,{t) 


P>{t) 


2 

4(<;>(r))2 


(23) 


where  <t>{r)  is  defined  in  Eq.  56. 

With  the  sine  wave  approximation,  the  duration  of  exceedance  of  a  threshold  can 
be  used  to  determine  the  amplitude  of  the  sliding  estimator  function,  provided  that  the 
angular  frequency  is  known.  We  assume  that  the  angular  frequency  can  be  adequately 
approximated  by  the  frequency  wj  corresponding  to  the  zero-crossing  period  Tj,  as 
illustrated  in  Fig.  8.  The  accuracy  of  this  iissumption  will  be  briefly  discussed  below. 
We  have 
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For  a  sinusoidal  exceedance  of  the  threshold  a  of  duration  t  with  frequency  u-m,  the 
amplitude  5,o  is  given  by 

a\  2{  1  —  cosojit) 

S,o  =  ^  (25) 

sinu;2T 

As  before,  the  threshold  for  the  starboard  sliding  estimator  function  is  given  by 

Putting  these  results  together,  we  have  the  following  numerical  procedure  for  esti¬ 
mating  the  probability  of  occurence  of  slides  of  distance  greater  than  an  arbitrary  length 
/.  We  determine  the  amplitude  of  the  sliding  estimator  function  5io  which  results  in 
a  slide  of  distance  I  by  solving  Eq.  22  iteratively.  The  duration  t'  resulting  in  this 
amplitude  is  obtained  by  solving  Eq.  25 


The  fraction  of  threshold  crossings  that  result  in  slides  uf  distance  greater  than  /  is 
given  by 

/“p,(rMr  =  exp(-j^(r‘)’)  (27) 

The  number  of  slides  per  second  of  distance  greater  than  /  is  obtained  by  multiplying 
Eq.  27  by  A/*,  which  is  given  in  Eq.  12. 

As  an  example  of  an  application  of  the  method,  we  will  revisit  the  chair-sliding 
experiment  considered  above  and  compute  the  expected  number  of  slides  of  distance 
greater  than  1.8  metres  (6  feet)  in  a  twenty  minute  period.  No  estimates  of  the  dynamic 
coefficient  of  friction  were  available,  but  for  the  purposes  of  the  exercise,  ^  was  assumed 
to  be  0.15.  As  before,  the  sliding  threshold  was  at  —  <t>m9  =  0.131^.  For  this 
experiment,  U2  =  0.679,  and  it  was  found  that  a  sliding  estimator  amplitude  of  0.149  ^’s 
resulted  in  a  sliding  distance  of  1.8  metres.  The  corresponding  duration  of  exceedance 
from  Eq.  26  was  1.47  seconds.  From  Eq.  27,  the  fraction  of  slides  of  distance  greater 
than  1.8  metre.s  was  found  to  be  0.54,  so  that  the  expected  number  of  slides  of  distance 
greater  than  1.8  metres  in  a  twenty  minute  period  was  24.5  x  0.54  =  13.2. 

In  order  to  get  some  idea  of  the  accuracy  of  using  the  zero-crossing  frequency  u>2  in 
determining  th  •  amplitudes  corresponding  to  a  threshold  exceedance  of  given  duration, 
the  threshold  ci  ossings  of  the  chair  sliding  experiment  were  analyzed.  For  each  crossing, 
the  duration  and  amplitu(  ?  of  the  sliding  estimator  function  were  measured,  and  the 
corresponding  Frequency  wj  was  determined  from  Eq.  26.  The  resulting  estimates  of 
u)i  had  a  mean  value  of  0.61  and  a  standard  deviation  of  0.12.  This  compares  with  an 
UJ2  value  of  0.6F;  computed  from  Eq.  24.  This  suggests  that  Eq.  25  may  be  sufficiently 
accurate  for  priictical  purposes,  but  further  investigation  is  required. 
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PREDICTION  OF  SLIDING  INCIDENCE  WITH 
NONZERO  LONGITUDINAL  ACCELERATION 


In  this  section,  we  present  a  frequency-domain  method  for  predicting  the  incidence 
of  sliding  in  the  presence  of  combined  vertical,  lateral,  and  longitudinal  acceleration. 
We  first  present  an  accurate  method  which  is  restricted  to  longitudinal  accelerations 
that  are  small  compared  with  the  lateral  acceleration,  w’hich  is  the  case  of  most  interest 
on  monohull  vessels.  W'e  then  present  an  approximate  method  for  dealing  with  the  case 
of  comparable  longitudinal  and  lateral  accelerations. 

We  first  consider  the  case  of  zero  heel  angle.  A  slide  will  occur  in  some  direction 
whenever 


y/Kng  +  K,  >  +  ^3) 

(28) 

v/^in,  +  F,L-w.i?3>/4? 

(29) 

Under  the  assumption  that  Fung  and  Fut  are  jointly  Gaussian,  it  can  be  shown 
that  y/F,o„j  -f  Fjo,  follows  the  Rayleigh  distribution.  Unfortunately,  the  distribution 

of  yJK.,  +  F,l.  -  nDs  is  the  convolution  of  a  Rayleigh  distribution  with  a  Gaussian 
distribution,  and  is  not  tractable.  It  is  therefore  necessary  to  develop  approximate 
solutions. 

We  first  consider  the  case  in  which  the  longitudinal  acceleration  is  small  compared 
with  the  lateral  acceleration.  Suppose  that  at  a  given  instant  in  time  the  resultant  force 
parallel  to  the  deck  (the  vector  sum  of  Fj,*  and  Ftong)  is  in  the  direction  6  shown  in 
Fig.  9.  The  resultant  force  will  result  in  a  slide  in  this  direction  provided  that 

Ftong  cos  6  +  Fiot  sin  6  -  /4D3  >  fj,g  (30) 

Under  the  assumption  that  Ftong  is  small  (but  not  necessarily  negligible)  compared  with 
Fia,,  most  of  the  sliding  incidents  will  take  place  within  some  range  about  S  =  90®  and 
6  =  270®,  and  slides  near  6  =  0®  and  6  =  180®  will  be  very  unlikely. 

We  will  estimate  the  number  of  slides  to  port  by  determining  the  incidence  of  up- 
crossings  of  by  the  resultant  of  Ftong  and  Sp  =  Ft^t  —  for  which  6  is  in 

the  neighbourhood  of  90®.  The  resultant  of  Ftong  and  Sp  is  in  the  direction  6*  = 
arctan(Fiaj  —  Ftong)  and  has  magnitude 

Ftong  cos  6*  -1-  Ftat  sin  6"  -  /4.D3  sin  6‘  (31) 

whereeis  the  direction  of  the  resultant  force  parallel  to  the  deck  is  given  by  6  = 
arctan(F/aj/F(onj)  and  the  condition  for  sliding  in  the  direction  S  is  given  by  Eq.  30. 
We  will  approximate  S  by  6*,  and  approximate  the  sliding  con  dition  by 
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Fieng  COS  6*  +  Ftat  sin  S'  ~  ^£>3  sin  S'  >  ^j,g 


(32) 


A  preliminary  assessment  of  the  accuracy  of  this  approximation  was  made  by  examining 
recorded  acceleration  time  histories  from  the  T-AGOS  monohull.  Equations  30  and 
32  were  both  evaluated  in  the  time  domain,  and  the  resulting  predictions  of  the  total 
number  of  sliding  incidents  were  compared.  From  these  results,  it  appears  that  Equation 
32  provides  a  good  approximation  to  Equation  30  provided  that  the  RMS  longitudinal 
acceleration  is  less  than  about  half  of  the  RMS  lateral  acceleration. 

As  an  aside,  we  note  that  the  method  presented  here  for  small  longitudinal  accelera¬ 
tion  would  apply  for  arbitrary  longitudinal  acceleration  if  the  vertical  acceleration  were 
zero.  For  monohulls,  it  is  difficult  to  conceive  of  a  situation  in  which  this  combination 
would  occur. 

We  assume  that  Fiong  and  Sp  are  jointly  Gaussian  with  zero  means,  standard  de¬ 
viations  Otong  and  <?p,  respectively,  and  correlation  coefficient  c  =  E{FiongSj,]  j atongf^p, 
where  E{FiongSp}  denotes  the  expected  value  of  FiongSp.  Under  these  assumptions 
their  joint  probability  density,  /{Fiong,  Sp)  is  given  by 


f{Eiong,Sp) 


_ 1 

27r<T((jnjOpN/T 


'/on  § 


(33) 


Although  Fiong  and  Sp  are  not  independent,  we  can  apply  a  coordinate  rotation  c  to 
obtain  independant  random  variables 


/  2  /  Ftong  cos  f  +  Sp  sin  c 

\w  )  y  —Fjons  sin  € -h  Sp  cose 


(34) 


provided  that  e  is  given  by 


tan  2c  = 


ong^p 

2 

^  long 


(33) 


W^e  will  estimate  the  number  of  slides  near  S  =  90®  as  the  number  of  nperossings  c>f 
the  tlireshold  n,g  by  the  resultant  of  z  and  w  in  the  half  plane  u»  >  0.  Slides  for  which 
S  is  not  in  the  neighbourhood  of  90®  will  be  rare  unless  Oig^g  is  a  significant  fraction  of 
H,g,  and  in  this  case  Op  will  be  sufficiently  large  to  cause  wide-spread  sliding  anyway. 

As  shown  in  Fig.  10,  the  number  of  slides  that  occur  at  a  distance  z  =  zq  from 
the  XL'  axis  can  be  estimated  as  the  sum  of  the  number  of  w  up-crossings  of  the  level 
\J{yL,gy  —  Zq  which  occur  in  the  small  interval  dzo  about  Zq  plus  the  number  of  c  up- 

crossings  of  the  level  ro  which  occur  in  the  small  interval  dwo  about  —  *o-  An 

estimate  of  the  total  number  of  slides  to  port,  Mp,  is  obtained  by  integrating  over 
Zq.  Denoting  the  ith  spectral  moments  of  z  and  w  by  m,  ^  and  respectively,  and 
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making  use  of  the  independence  of  these  two  random  variables,  we  obtain 


°‘d:o  + 


(36) 


-  4 

+ 


where 


b  = 


\  ^0,r  ^0,u>  / 


(37) 


and  ‘erf’  denotes  the  error  function,  which  is  defined  in  Eq.  58  of  Appendix  A.  The 
integral  over  y  in  Eq.  36  appears  to  be  improper;  however,  making  the  substitution 
^  I  ^  v;e  obtain 


Jo  Jo 


(38) 


In  general,  this  last  integral  can  be  evaluated  numerically;  however,  in  the  case  that  b  is 
large,  most  of  the  contribution  to  the  integral  over  y  in  Eq.  36  will  come  from  the  points 
near  the  origin.  In  this  case,  the  integral  can  be  estimated  by  expanding  l/\/l  — 
about  y  =  0.  We  obtain 


i: 


y^y.  1  .  -L  _  4.  J- 


\/i  -  y* 


2b  46* 


46  462 


(39) 


The  spectral  moments  of  z  and  w  can  be  expressed  in  terms  of  the  moments  and 
cross-spectral  moments  of  Fiong  and  Sp.  Let  (Fionp.o/Co)  and  (5p, o/Co)  denote  the  (coirr 
plex)  transfer  functions  of  Fiong  and  5p,  respectively,  and  5^(0;^)  denote  the  wave 
trum  as  a  function  oi*  encounter  frequency  w*. 


mo 


-=r( 


F long.O  ,  *^p,0 

— COS  c  +  sin  e 

Co  Co 


cos  €  +  ^- sin  e  1  S((u)*)du)«  (40) 
Co  Co 


where  the  bars  indicate  complex  conjugates.  We  have 

^o.z  =  cos^  €  + 2E{F/on3Sp}  sine  cose  +  o-p  sin^  €  (41) 

=  cr^ong  cos^  e  +  2cffiarig<yp  sin  e  cos  e  +  sir^  e 

Making  use  of  Eq.  35,  we  obtain 
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(42) 


mo.:  -  ^{(f^ong  +  “^p)  “  9\/(‘^p  ~  ^^/onp)^  +  (2ca,on5ap )* 

Similarly, 

mo.u,  =  2(<^L5  +  <7p)  +  +  (2c<T,on,<7p)*  (43) 

The  second  moments  can  be  determined  in  an  analogous  fcishion.  If  we  denote  by  I  the 
integral 


^long.Q  ^p.O  .  ^ [on^.O  *^p,0 

Co  Co  Co  Co  , 


S(:(We)du;, 


we  obtain 


^  _  I  _  V  (”^2,P  m2.iong){<^l  <^}ong) 

m2,:  =  -[rn2,tong  +  JTlZ.p) - 1 - ^  ^ 

^  -  (^Jor^gY  +  i2ca,ongCTpy 

I C  ^ tong^p 


\/K-^?onp)^  +  (2c«T„„,CTp)^ 


and 


^  -  <^Lp)*  +  {2cai,ng(Jj,y 

IC(JlQf^g(J'p 


-(^Ing^  +  (2c<7,ong<ypy 


(44) 


(45) 


(46) 


The  incidence  of  sliding  in  the  starboard  direction  can  be  determined  in  a  similar 
fashion.  In  the  presence  of  a  mean  heel  angle  (prn  the  threshold  (j,g  should  be  replaced 
by  +  4>m9  ^  ~  4>mg,  for  port  and  starboard  sliding,  respectively. 

To  gain  some  insight  into  the  increased  sliding  incidence  which  resuits  from  the 
presence  of  nonzero  longitudinal  acceleration,  we  will  compare  the  port  sliding  incidence 
predicted  from  Eq.  36  with  the  sliding  incidence  from  Eq.  12  for  a  numerical  example. 
We  assume  Ciong  =  <^p/^  and  cTp/S  =  0.5p,<7,  c  =  0.5,  and  <f>m  =  0.  For  simplicity,  we 
take 


{l/2Tr)^Tn2,p/mo,p  =  (l/2Tr)^m2,^lmo,^  =  (l/27r)^m2.,/mo,z  (47) 

For  this  example,  we  find  that  the  presence  of  longitudinal  acceleration  increases  the 
sliding  incidence  by  about  14%. 

For  longitudinal  acceleration  larger  than  about  half  the  lateral  acceleration,  it  is  no 
longer  reasonable  to  assume  that  all  of  the  slides  will  occur  near  the  port  and  starboard 
directions.  An  estimate  of  the  number  of  sliding  incidents  in  the  case  of  comparable 
lateral  and  longitudinal  accelerations  can  be  obtained  as  follows.  The  coefficient  of 
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friction  ^  is  the  same  at  all  angles  of  azimuth;  hence,  a  polar  plot  /i,  as  a  function  of 
8  forms  a  circle.  We  will  (somewhat  crudely)  approximate  this  circle  by  a  square  of 
equal  area,  as  shown  in  Fig.  11.  The  coefficient  of  friction  now  becomes  a  function  of 
8,  Suppose  that  at  a  given  instant,  the  resultant  force  parallel  to  the  deck  is  in  the 
direction  8  as  shown  in  Fig.  12.  The  condition  for  sliding  in  the  direction  8  which  was 


given  in  Eq.  30  now  becomes 

Fiong  cos  8  +  Fiat  sin  8  -  ^{8)b:^  >  ^{8)g  (4^.^) 

Note  however,  that  in  order  for  a  slide  to  occur  in  some  direction  8  in  the  quadrant 
within  45®  of  the  longitudinal  direction,  it  is  necessary  and  sufficient  that 

Flong  -  /4(0)D3  >  /4(0)y  (49) 

To  see  this,  w^e  note  that  Eq.  48  is  equivalent  to 

(50) 

Multiplying  this  result  by  cos^  =  + n,  we  obtain 

Flong  fls{8)  COS  8b3  >  fJ^{8)  COS  8  g  (51) 

w'hich  is  equivalent  to  Equation  49. 

If  we  define  the  forward  sliding  estimator  function  as 

Sf  =  Flong  -  /4(0)D3  (52) 

then  the  forward  slides  can  be  estimated  by  counting  the  upcrossings  of  5/  of  the 
threshold  ^{0)g,  Similarly,  we  define  the  aft  sliding  estimator  function  via 

Sq  =  -Flong  -  fjk{0)b3  (53) 


The  total  number  of  slides  can  be  approximated  by  summing  the  upcrossings  of  the 
sliding  estimator  functions  in  the  four  primary  directions.  For  the  equal  area  square 
approximation  to  the  circle,  we  have  /i»(0)  =  12  =  0.886;ig. 

The  reader  is  cautioned  that  in  the  case  of  small  longitudinal  acceleration,  this 
approximation  results  in  an  overestimate  of  the  number  of  slides.  For  the  example  con¬ 
sidered  above  in  which  criong  =  <^p/3  and  ^p/3  ==  the  number  of  sliding  incidents 

predicted  to  occur  to  port  with  the  reduced  threshold  0.886/i,  is  54%  larger  than  the 
sliding  incidence  with  the  original  threshold.  The  more  careful  analysis  presented  above 
showed  that  the  presence  of  longitudinal  acceleration  produced  an  increaise  of  only  14%. 
The  estimates  obtained  w'ith  the  square  approximation  should  improve  in  the  caise  of 
comparable  longitudinal  and  lateral  accelerations. 
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CONCLUDING  REMARKS 


A  frequency‘doma*n  method  of  predicting  the  incidence  of  personnel  or  equipment 
sliding  has  been  presented.  The  method  includes  the  (linearized)  forces  due  to  roll, 
pitch,  longitudinal,  lateral  and  vertical  accelerations,  and  also  the  effects  of  non-zero 
mean  heel.  Two  different  solutions  were  developed  according  to  the  relative  size  of  the 
longitudinal  acceleration.  The  method  developed  for  small  longitudinal  acceleration  is 
accurate,  while  the  method  developed  for  comparable  longitudinal  and  lateral  accelera¬ 
tions  is  more  approximate,  but  should  still  produce  estimates  of  sufficient  accuracy  for 
practical  purposes.  The  predictions  of  the  method  were  compared  with  the  results  of  a 
chair-sliding  experiment  conducted  at  sea,  and  good  correspondence  v/bls  obtained. 

A  simple  approximation  for  the  duration  of  threshold  exceedance  probability  density 
w^as  described,  and  used  to  develop  a  method  for  predicting  the  incidence  of  slides 
greater  than  a  certain  severity.  The  method  applies  for  arbitrary  lateral  and  vertical 
acceleration,  but  only  to  the  case  in  which  the  longitudinal  acceleration  is  negligible. 
An  example  of  the  application  of  the  method  was  presented. 

It  has  been  proposed  that  in  future,  seakeeping  criteria  should  be  presented  in  terms 
of  the  tolerable  incidence  of  degrading  events  such  as  helicopter  sliding.  Criteria  in  this 
form  are  simple  and  are  equally  applicable  to  all  vessel  types. 
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i.  Starboard  Sliding  Estimator  Functions  and  Thresholds  Correspond 
ing  to  Different  Coefficients  of  Friction 
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Fig.  7.  Sliding  Estimator  Function  in  Regular  Waves 
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Fig.  9.  Definition  of  Sliding  Direction  6 


Fig.  10.  Illustration  of  Sliding  Calculation  for  Small  Fung 
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Fig.  11.  Ekiual-Area  Square  Approximation  to  the  Circle^ 
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Fig.  12.  Sliding  Calculation  in  the  Case  of  Longitudinal  and  Lateral  Accel¬ 
erations  of  Comparable  Magnitude 


APPENDIX  A 

A  SIMPLE  APPROXIMATION  FOR  THE  DURATION 
OF  THRESHOLD  EXCEEDANCE  PROBABILITY  DENSITY 

by 

Ross  Graham  and  John  F.  Dalzell 


INTRODUCTION 

In  seakeeping  research,  there  are  frequent  occasions  when  the  number  of  up-crossings 
of  a  threshold  by  a  zero-mean,  Gaussian  stochastic  process  must  be  determined.  Stan¬ 
dard  examples  include  estimating  the  probability  of  keel  emergence  and  the  probability 
of  deck  wetness.  The  solution  to  this  problem  is  well-known,  and  was  first  developed 
by  Rice*  (see  also  Price  and  Bishop*). 

For  certain  applications,  information  on  the  duration  of  time  that  the  process  spends 
above  the  threshold  is  also  required.  Complete  information  on  the  duration  statistics 
would  be  provided  by  the  duration  of  exceedance  probability  density,  but  unfortunately, 
this  density  can  only  be  estimated  numerically,  and  the  available  approximations  are 
both  difficult  to  implement,  and  suffer  from  serious  problems. 

The  first  of  the  available  approximations  originates  with  the  theory  of  Rice",  as 
refined  by  Kuznetsov,  Stratanovicb  and  Tikhonov"'**.  (The  theory  and  approximations 
are  also  partially  described  by  Price  and  Bishop".)  The  “exact”  theory  deals  with  the 
conditional  probability  density  for  the  duration  of  exceedance  of  an  arbitrary  threshold, 
with  no  pre-conditions  about  the  nature  of  the  process  or  the  statistical  independence 
of  the  times  of  threshold  crossings.  The  conditioning  of  the  density  is  upon  the  value 
of  the  threshold,  with  a  duration  event  defined  as  the  time  interval  between  successive 
up-  and  down-crossings  of  the  threshold.  This  “exact”  theory  takes  the  form  of  an 
infinite  functional  series  which  was  impossible  to  evaluate  in  the  IQSO’s — and  probably 
still  is.  Kence,  practical  use  has  involved  approximations  and  assumptions  of  one  sort 
or  another. 

The  boldest  simplifying  assumption  is  that  c//  the  times  of  threshold  crossings  are 
statistically  independent.  With  this  assumption,  the  durations  become  a  Poisson  pro¬ 
cess,  a  model  whicl  appears  relevant  for  very  high  thresholds,  and  which  has  been  noted 
in  the  seakeeping  literature*^. 

The  usual  approach  for  arbitrary  thresholds  is  due  to  Rice**  and  Tikhonov**,  and  will 
be  called  “Approximation  1”  herein.  This  approximation  involves  taking  only  the  first 
term  of  the  infinite  series  of  the  exact  theory  as  an  approximation  for  arbitrary  threshold 
levels  but  “small”  durations.  Unfortunately,  this  approach  yields  a  result  which  is  not 
exactly  that  desired.  In  effect,  the  approximate  density  involves  the  probabilities  that, 
given  a  threshold  up-crossing  at  some  time,  there  is  a  down-crossing  threshold  in  the 
neighborhood  of  a  later  time  t.  regardless  of  what  happened  between  time  fo  and  to  +  t. 
As  a  result,  Approximation  1  appears  to  yield  good  approximations  to  the  true  density 
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for  short  dur&tions,  but  does  not  integrate  to  unity,  and  in  fact  tlie  integral  may  not 
even  be  finite 

An  essentiilly  more  modern  treatment  of  the  problem  is  due  to  Rainal*^  and  Mi¬ 
maki'®’’^.  In  these  works,  the  estimation  of  the  conditional  density  was  formulated 
in  a  different  way,  and  by  making  the  assumption  that  a  given  individual  interval 
between  thres  lold  crossings  is  statistically  independent  of  nearly  all  the  possible  sums 
of  preceding  iiitervals,  the  problem  was  reduced  to  the  solution  of  a  Volterra  integral 
equation  involving  the  numerical  results  from  Approximation  1.  described  above.  On 
the  whole.  Approximation  2,  though  more  computer  intensive,  yields  more  satisfying 
results  than  Approximation  1.  However,  the  approximation  has  disquieting  faults  of  its 
own — principjlly  that  the  computed  density  may  become  negative  for  large  durations. 

Thus,  both  of  the  available  approximations  for  arbitrary  threshold  levels  are  useless 
when  moment  ?  of  the  density  are  required.  The  purpose  of  this  appendix  is  to  examine  a 
high-threshold  approximation  which  was  developed  in  Reference  11,  but  which  appears 
to  have  gone  unnoticed.  The  approximation,  which  is  described  in  the  next  section,  has 
a  simple  analytic  form  and  satisfies  the  basic  requirements  of  a  probability  density.  The 
question  of  how  high  the  threshold  must  be  before  a  good  approximation  is  obtained  is 
investigated  ir  a  later  section  using  the  results  of  a  time  domain  simulation  of  a  typical 
long-crested  wave  system. 


APPROXIM  ATION  FOR  THE  DURATION  OF  THRESHOLD 
EXC.EEDANCE  PROBABILITY  DENSITY 

Let  St{t)  be  a  zero-mean,  Gaussian,  stochastic  process  (not  necessarily  narrow- 
banded,.  In  the  main  text,  Sg{t)  is  assumed  to  be  the  starboard  sliding  estimator 
function,  but  is  an  arbitrary  function  in  this  appendix.  Let  a,  denote  the  standard 
deviation  of  5,(<),  and  suppose  that  the  threshold  level  is  at  a.  We  define  the  nondi- 
mensional  threshold  by  r  =  a/<7*. 

We  will  denote  by  r  the  duration  of  a  given  threshold  exceedance,  and  the  duration 
probability  demsity  by  Rice’*  demonstrates  that  for  high  thresholds,  p,{t)  tends 

towards  the  Rayleigh  probability  density,  i.e. 

p,(r)  -»  (54) 


Rice  (see  also  Vanmarcke’*)  shows  that  the  average  duration  f  of  a  threshold  exceedance 
is  given  by 


w’here  M,  is  given  in  Eq.  12,  and 


(55) 


(56) 
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Equation  55  is  not  restricted  to  high  threshhold  levels.  The  function  d(r)  is  related  to 
the  ervor  function  erf(i/)  via 


For  the  Rayleigh  probability  density,  we  have 

f  =  =  /Z  (59, 


so  that 


A  = 


•kM] 


2f^  l[4>{r)Y 

This  results  in  the  following  approximation*’  for  p*(r) 


(60) 


(61) 


for  high  enough  thresholds. 


TIME  DOMAIN  SIMULATION 

In  order  to  get  some  feeling  for  how  high  a  threshold  must  be  in  order  to  obtain 
reasonable  results  with  the  Rayleigh  approximation  of  Eq.  61,  a  typical  long-crested 
wave  system  was  modeled  as  a  zero-mean  Gaussian  process,  simulated  in  the  time 
domain,  and  the  resulting  “data”  used  to  derive  some  duration  statistics  for  comparison 
with  the  approximation. 

Though  the  results  were  later  nondimensionalized,  the  variance  spectrum  of  the 
long-crested  wave  system  was  chosen  to  be  of  the  Bretschneider  type  with  modal  pe¬ 
riod,  Tm,  of  9  seconds  and  significant  height  of  4  units.  The  mechanics  of  the  time 
domain  simulation  involved  the  “fast  convolution”  method  in  which:  pseudo  random 
bandlimited  white  noise  is  Fast  Fourier  transformed;  the  resulting  transform  is  mod¬ 
ified  by  a  nonrealizable  frequency  domain  filter  corresponding  to  the  desired  variance 
spectrum;  and  an  inverse  Fast  Fourier  Transform  then  in  turn  yields  the  simulated 
process. 

The  time  series  simulations  of  realizations  of  the  selected  wave  spectrum  were  gener¬ 
ated  in  “handy-sized”  realizations  of  8K  point  time  series,  at  delta-time  of  0.1  seconds. 
Thirty  independent  realizations  of  the  process  were  simulated  by  entering  the  computer 
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noise  generator  at  widely  separated  points  in  its  sequence.  Considering  all  of  the  sam¬ 
ples,  the  simulation  was  equivalent  to  a  6.7  hour  observation  of  the  wave  process  defined 
by  the  9  second  modal  period  spectrum. 

A  threshold  crossing  analysis  was  performed  on  the  simulated  data  for  six  nondi- 
mensional  thresholds  corresponding  to  r  —  0,  0.5,  1.0,  1.5,  2.0  and  2.5.  The  mean  time 
interval  between  threshold  crossings  wtis  first  computed,  and  found  to  agree  well  with 
the  corresponding  theoretical  estimates  (1/A/j),  Eq.  12.  The  computed  time  differences 
between  each  up-  and  down-crossing  were  combined  into  an  ordered  sample  for  further 
analysis,  and  each  duration,  r,  was  nondimensionalized  by  the  wave  modal  period.  T^- 

For  purposes  of  comparisons  with  theory,  an  estimate  of  the  mean  probability  density 
over  some  small  interval  of  duration  was  required.  A  uniform  nondimensional  duration 
class  interval,  Alr/Tm},  was  first  assumed  and  the  number  of  simulated  exceedance 
durations  falling  into  each  was  tallied.  The  resulting  tallies  were  divided  by  the  total 
number  of  exceedance  durations  in  the  sample  to  form  an  estimate  of  the  probability, 
[p,(T/Tm)  A{r/r„,}],  and  then  this  result  was  divided  by  A{T/Tm}  to  form  an  estimate 
of  the  mean  probability  density  over  the  interval.  For  graphical  purposes  these  estimates 
were  associated  with  the  mid-point  of  the  nondimensional  class  interval. 

RESULTS  AND  DISCUSSION 

The  results  from  the  simulation  are  compared  with  Eq.  61  in  Figs.  13  through  18  in 
terms  of  the  nondimensional  duration  time  r/Tm.  Computations  of  duration  densities 
according  to  Approximations  1  and  2  as  described  in  the  introduction  are  shown  for 
comparison. 

As  expected,  the  comparison  is  poor  for  a  0  threshold  (Fig.  13),  and  at  most  fair  for 
r  =  0.5  (Fig.  14),  but  the  results  are  surprisingly  good  for  the  higher  threshold  levels. 
The  comparison  shown  in  Fig.  15  for  r  =  1.0  is  already  quite  acceptable,  and  Figs.  16, 
17  and  18  show'  that  good  agreement  is  obtained  a.  higher  thresholds. 

A  closer  examination  of  Figs.  13  to  18  reveals  that  the  average  duration  f  of  a 
threshold  exceedance  is  slightly  underestimated  in  all  cases.  In  theory,  Eq.  55  should 
hold  exeictly  for  all  threshold  levels.  The  discrepancy  appears  to  arise  from  differences 
in  the  theoretical  estimate  of  rrxj  in  Eq.  12,  and  the  value  of  mj  realised  in  the  simu¬ 
lated  time  history.  The  theoretical  estimate  was  obtained  by  direct  integration  of  the 
analytical  expression  for  the  Bretschneider  spectrum,  and  is  slightly  higher  than  the 
value  for  the  simulation.  This  discussion  serves  to  highlight  that  accurate  estimates  of 
the  threshold  exceedance  probability  density  depend  strongly  on  accurate  estimates  of 
mj. 

It  is  concluded  that,  at  least  for  the  sample  wave  elevation  process  considered  here, 
Eq.  61  provides  a  reasonable  engineering  approximation  to  the  threshold  exceedance 
duration  probability  density  function  for  nondimensional  thresholds  greater  than  or 
equal  to  one. 

It  should  be  noted  that  the  wave  elevation  process  is  nar'c.w-banded  The  applica- 


28 


bility  of  the  Rayleigh  approximation  for  the  exceedance  duration  probability  density  of 
wide-band  processes  requires  further  investigation. 

CONCLUDING  REMARKS 

A  simple  approximation  due  to  Rice"  for  the  duration  of  threshold  exceedance  prob¬ 
ability  density  at  high  threshold  levels  was  presented,  and  its  applicability  to  the  wave 
elevation  process  was  investigated,  by  comparison  with  the  results  of  a  time-domain 
simulation.  Unlike  other  more  complicated  approximations  in  use,  this  approximation 
has  a  simple  analytic  form,  namely,  the  Rayleigh  probability  density. 

It  was  found  that  the  simple  approximation  compared  well  with  the  time-domain 
simulation  results  for  nondimensional  thresholds  greater  than  or  equal  to  one. 
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P,(t/T 


r  s  0.0,  3603  Simulated  Events 


Fig.  13.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Nondimensional  Threshold  r  =  0.0 
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P.(t/T 


r  B  0.5,  3200  Simulated  Events 


t/T„ 


Fig.  14.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Nondimensional  Threshold  r  =  0.5 
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P,(t/T 


r  >=  1.0,  2230  Simulated  Events 


Fig.  15.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Nondimensional  Threshold  r  =  1.0 
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P,(t/T 


r  —  1.5,  1211  Simulated  Events 


Fig.  16.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Noodimensional  Threshold  r  =  1.5 
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P,(t/T 


r  >=  8.0,  506  Simulated  Events 


Fig.  17.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Nondimensional  Threshold  r  =  2.0 
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P,(t/T 


r  «  2.5,  150  Simulated  Events 


Fig.  18.  Comparison  of  Simulated  Threshold  Exceedance  Probability 
Density  Function  with  Approximations  1  and  2  and  the  Rayleigh 
Approximation  for  a  Nondimensional  Threshold  r  =  2.5 
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